Geodesic deviation in pp-wave spacetimes of quadratic curvature 

gravity 

Edgard C. de Rey NetcQ 

Instituto Nacional de Pesquisas Espaciais - Divisdo de Astrofwica 
Av. dos Astronautas 1758, 
Sao Jose dos Campos, 
12227-010 SP, Brazil 
(Dated: September 05, 2003) 

Abstract 

We write the equation of geodesic deviations in the spacetime of pp-waves in terms of the 
Newman-Penrose scalars and apply it to study gravitational waves in quadratic curvature gravity. 
We show that quadratic curvature gravity pp-waves can have a transverse helicity-0 polarization 
mode and two transverse helicity-2 general relativity-like wave polarizations. A concrete example is 
given in which we analyze the wave polarizations of an exact impulsive gravitational wave solution 
to quadratic curvature gravity. 
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I. INTRODUCTION 



In the present days, we assist a increasing worldwide activity to detect gravitational 
radiation. These efforts will certainly lead to a birth of the gravitational wave astronomy in 
the beginning of the present century. Consequently, there is a growing theoretical interest 
in the description of gravitational waves emission and detection. Besides the information 
on astrophysical process that can be gained from gravitational wave observations, the direct 
detection of gravitational waves will be also of great importance to the understanding of the 
nature of gravity. Within this context, although the common belief that the linearized 
theory is sufficient to describe gravitational waves, exact solutions of gravity equations 
which represents gravitational waves also deserves to be investigated. One reason to this 
is that the linearization of gravity equations can hide important features contained in the 
exact nonlinear systems which could be useful to the understanding the global content of 
the gravitational theories involved. For instance, exact radiative solutions to nonlinear 
curvature gravities are easily obtained for pp-wave metrics while, the contribution of the 
nonlinear curvature terms to linearized quadratic gravity waves are obscured by the linear 
approximation Moreover, as pointed out in j3|, the nonlinear content of Einstein's 

gravity cannot be neglected in many physically interesting cases such as the gravitational 
emission in the coalescence of binary systems. Also, in ^, exact solutions to general 
relativity which represents spin-1 gravitational waves are found. 

In the present work we study the geodesic deviations in pp-wave spacetimes in the 
framework of quadratic curvature gravity. The plane fronted gravitational waves with 
parallel rays, pp-waves, are spacetimes which admits a covariantly constant null vector 
field. These spacetimes represents plane gravitational waves which propagates with the 
fundamental velocity c. They constitute a subclass of the general Kundt class of exact 
plane gravitational waves [^. Here, we are interested in obtaining the effect of quadratic 
gravity waves on geodetic test particles. Then, we study the geodesic deviations in pp-wave 
spacetimes which we know to be simple exact solutions to quadratic curvature gravity 0|. 
The existence of exact pp-wave solutions to quadratic gravity leads to the question of 
the existence of solutions for the general Kundt class of exact plane waves and also for 
the Robinson- Trautman class of exact spherical gravitational waves in quadratic gravity. 
However, the investigations of this subject are beyond the scope of the present work. The 
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choice of the pp-wave spacetime metric also excludes the non-null wave like solutions of 
quadratic gravity, such as those found in [3]. 

The equation of geodesic deviations gives the relative accelerations between free test 
particles falling in a gravitational field and is a cornerstone to the understanding of the 
physical effects of the gravitation j^, being the basis of almost all prospects for the 
gravitational wave detection. Geodesies and geodesic deviations in spacetimes of impulsive 
pp-waves of general relativity are rigorously studied in by using the concepts of the 
Colombeau algebras to handle the nonlinear products of distributions. We do not follow 

n 

here the same approach developed in 9'], rather than, we write the equation of geodesic 
deviations in a orthonormal tetrad basis by projecting the components of the Riemann 
tensor on this basis. Then, we obtain the relative accelerations of nearby test particles as 
function of the Newman-Penrose (NP) quantities. Another serie of recent papers which deals 
with geodesic deviations in spacetimes of general relativity and, in particular with the field 
of pp- waves deserves to be cited here [l^ • However, we stress that our approach to the issue 
of geodesic deviations in the field of pp-waves also differ from that carried out in the above 
references. 

The plane of the paper is the following. In the section |n] we write a NP null tetrad in 
the spacetime of generic pp-waves and the nonvanishing NP quantities. In the section IIIII 
we write the geodesic deviation equation in terms of the NP quantities. Then, we reduce 
to the case of pp-waves and obtain the wave polarizations which can arise for these metrics. 
In the section IIV| by considering pp- waves which are solution of non-linear Lagrangian field 
equations, we obtain the geodesic deviations as such as, the wave polarization that can be 
found in quadratic curvature gravity. We also analyze a concrete example given by the 
gravitational wave solution to quadratic curvature gravity obtained in ^|. 



II. NEWMAN-PENROSE FORMALISM FOR GENERIC pp- WAVES 

In this section we extend the method employed by Hayashi and Samura in to construct 
a null tetrad frame for cylindrically symmetric gravitational shock waves, to the case of waves 
with arbitrary amplitude and arbitrary dependence on angular coordinates. Then, we write 
the non-necessarily null NP quantities for a generic pp-wave. 



3 



We start with a pp-wave described by the hne element spacetime {h = c = 1): 

ds^ = —dudv + H{x'', u)du^ + gkidx^dx\ (1) 

where u = t — z,v = t + z and are generahzed coordinates in the 2- dimensional space 
perpendicular to the propagation direction, the so called transverse space. The Latin indices 
run over the transverse space dimensions k = 1,2. The line element (PQ) clearly represents 
a plane wave propagating in the z direction with the fundamental velocity. Let us write 
H = H{x±, u) where x± is a point in the transverse space. The particular form of H{x±, u) 
will depend on the source term as such as on the underling theoretical model. 

The non-necessarily null Christoffel symbols for the spacetime defined by are 



dH 

rL = -||(^^,«) , ri = -i,^'|^(x^,n), (2) 

where f are the transverse space Christoffel symbols. In the spacetime (^J we have 

R = ; Rf,^R^'' = and R^^apR^"''"^ = 0, (3) 

and the nonvanishing components of the Ricci tensor are 

Rki = Rki and Run = -^'^]_H{x^,u), (4) 

where V\ is the Laplacian in the transverse space. In fiat transverse space we have Rki = 0. 
The geodesic equations are 

ii = 0, (5a) 

■■k kl^'^ -2 T^k •I'm /"c \ 

where the dot means derivative with respect to the affine parameter. 

Now, we consider flat transverse space and define cylindrical coordinates on it by writing 
x^ = pcos(0), = psin(0). This choice makes apparent the radial and angular aspects of 
the pp-wave metric. The line element (0) takes the form 



ds^ = -dudv + H{p, 0, u)du^ + dp^ + p^dcfr. (6) 



The geodesic equations becomes 

u = 0, 



(7a) 



.. OH dH.. ^dH . 

v = -u^ + 2-pu + 2-^u, (7b) 



I OH 2 



e - -(j)p, (7c) 



2p2 d(f) p 

Let us write a point in the spacetime with coordinates {u, v, p, 0) as x. To construct a 
null tetrad basis we need to know the first derivative of the x with respect to the affine 
parameter. This is obtained by integrating the geodesic equations once with respect to the 
affine parameter and holding the transverse coordinates fixed. 

The first of equations ((Tj) implies that u = aos + bo. Without loss of generality, we choose 
6o = and Oq = 1 so that u can be taken as the affine parameter itself. Thus, 

u = l, (8a) 



v = 2H- I —du, (8b) 

/no 



1 r dH 



2p 



du, (8c) 
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where we neglect the inessential constants. 

Now, we define the vector / by the vector whose contravariant components are 
/'^ = (-it, p, 0), /X = 0, 1, 2, 3. By noting that 

p2 + p^' = H- r ^du, (9) 
Juo du 

(take the total derivative with respect to u of both sides of Q) one can show that l^l^ = 0. 

The equations (jHl) defines a null tangent vector to the spacetime ©. We can build a 
Newman-Penrose tetrad by taking the vector I above, a null real vector k orthogonal to I 



and two null complex conjugated vectors m and m required to satisfy the orthogonality 
conditions: 

lm = lni=krn = krri = 0, (10) 

the null conditions: 

ll = kk = rnrn = ffirri = (11) 

and the normalization conditions: 

k ■ I = —1 and m ■ m = 1. (12) 

The null vector k must be proportional to S^a; P| . Taking into account the normalization 
conditions (fT^ we obtain 

A;^ = (0,2,0,0) ; A;^ = (-1, 0, 0, 0). (13) 
For the vector m we take 

m^ = (0,«,/?,7), (14) 

where a, P, and 7 may be complex. The conditions (fTUj) . (fTTjl and (fT^ determines a, /3 and 
7. Thus, 

« = ^(p + ^P0) , /? = ^ , 7 = ^. (15) 

By following the standard notation, we can write the null tetrad basis {e(a)} = 
{l,k,m,m}, 6(1) = l^^ , e(,) = k^^ , 6(3) = m/^l, , 6(4) = m-^, and the dual 
basis {e(a)}, e(i) = l^d^f" , 6(2) = fc^.ci^'' , 6(3) = m^d^'' , 6(4) = m^di^, where 

F = {u,v,p,4'), (16a) 



2/ 9m 2' 



F = (0,2,0,0), (17a) 



A;^ = (-1,0,0,0), (17b) 
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= (^0, 2(p + ip^), 1, , (18a) 

= |^-^(p + ip0),O,l,ip^ . (18b) 

Note that the components are the momenta canonically conjugated to the coordinates 
(m,w,P,0). 

Using the equations ()16|) - (fTH|) and the definitions in the Appendix A, one can verify that 
the only non-necessarily null NP quantities are: 

$22 = = -\vlH (19) 



and 

= C,.-,d>-nfPnf = i (^ + 5£ _ H,^ + 1 (^ - :^) , (20) 

where Cf^^jS is the Weyl tensor and the partial derivatives are abbreviated by a comma. We 
can immediately see that if if is a harmonic function of the transverse coordinates, $22 = 0. 
Moreover, if H is cylindrically symmetric, 53m\E'4 = 0. 

III. GEODESIC DEVIATION OF TEST PARTICLES 

In this section we write the components of the Riemann tensor in terms of the NP 
quantities in a local orthogonal basis. We follow the approach of , but we do not impose 
any field equation to the spacetime metric. Thus, we obtain a model independent description 
of geodesic deviation equation in pp -wave spacetimes. 

Let u = d'^/dr, u^u^ = —1 be the four- velocity of free a test particle in a spacetime of 
curvature described by the Riemann tensor R^v'-iS- The displacement vector X between two 
nearby particles must obey the geodesic deviation equation 

= -W,-,sv!'X"<u', (21) 

where r is the proper time of one of the particles. 

We define a local basis by the orthonormal tetrad {ea}, a = 0, 1, 2, 3, where = u and 
{e-}, i = 1,2,3 are orthogonal space-like unit four-vectors, such that ea ■ = g^^e^e^ = 
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r^-g = diag(— 1, 1, 1, 1). The frame components are X"- = e^X'^. The relative accelerations 
between two test particles in the local basis are defined by 



Then, we have 



■ -rr- (22) 



dr2 
and 



^^^^ R^,~,5u''u''X^u^ = (23) 



d^X' 



dr^ ' 
where X*=e^X'^ gives the distance between two test particles and 



R'mX^ (24) 



^'iOjo ~ ^fj-i^-yS^i^o^l^i (25) 



are the projections of the Riemann tensor components on the local basis {e^}. The local 
basis is related to a null tetrad basis {I, k, m,m} by ^: 

u=^{l + k) , eg = (fc - Z) , 

ei = ^(m: + m) , = i-^ (m - m) . (26) 

From the definition of the Weyl tensor we have 

-^iojo = Qojo + 2^^ijRd6 - ^ij) - Q^ijR- (27) 

By using ()26|) and ()27|) . we can write the components of Riemann tensor on the local basis 
in terms of the twelve NP quantities. The result is in the Appendix A, equations ()A4jl . For 
the pp-wave spacetime the only non- necessarily null NP quantities are $22 and \E'4. Thus, 
the equations ()A4|) reduces to: 

i?ioio = ^3?evl>4 + ^$22, (28a) 



^1626 = -T^^rn^i , (28b) 



^202o = -^3f^evI/4 + ^$22. (28c) 
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The equations (|2^ reads 



d^X^ 



-(A+ + Ao)X^ + AxX2, (29a) 



A^X' - {-A+ + A,)X\ (29b) 
2 0, (29c) 



df^ 

d^X^ 



dt 

where A+ = |3?e^'4, A^ = |5>m^4 and Ao = i$22- 

As can be immediately seen from ()29|1 . the generic pp-wave produces no variation in the 
longitudinal direction Cg. Under a rotation of the transverse plane by an angle i), m changes 
according to m' = e~^^m. Then, 

= cos(2^9) - sin(2'^) Ax , (30) 



a; = cos(2i9)Ax + sin(27?)A+, (31) 

K = Ao. (32) 

By using the above equations, it is easy to see that the equations (j^^ are invariant under 
rotations of ^9 = rnr, where n is an integer. As Ao is invariant under arbitrary rotations 
of the transverse plane we conclude that A^ and Ax are responsible to the "+" and "x" 
helicity-2 polarization modes and Ao to the helicity-0 mode. By a rotation of -(9 = n7r/4 
we have A'_^ = —A^ and A'y, = A+. Then, a general observer sees a superposition of two 
helicity-2 polarization modes shifted by 7r/4 and one helicity-0 polarization mode. From (fTUj) 
we can see that, if V±H = 0, there are no helicity-0 polarization modes. As can be seen 
from ()20|) . if the transverse space is cylindrically symmetric, A^ = and the wave is purely 
polarized. Note that, we have not made any assumptions concerning the gravitational 
field equations or the underling theoretical model. The results obtained until now follows 
directly from the structure of the pp-waves. 

We also notice that, the line element (or equivalently (jH))) is of Kerr-Schild form P|. 
This means that, the spacetime metric Qf^u can be written as 

9fM,y = ri^,y + Hp^,py, (33) 



where = —5^u- Then, the weak field regime is obtained simply by taken H <^ 1. Note 
that since the pp-waves are exact solutions to general relativity and even to more general 
g.av.a*„^.heoHeJ. the weak aeld regime described here are also exact solutions since, 
H are not imposed to satisfy the linearized gravity equations. 



IV. WAVE POLARIZATIONS IN QUADRATIC CURVATURE GRAVITY 

In this section, we apply the formalism developed in the two previous sections to analyze 
the relative acceleration of test particles in presence of pp-waves in quadratic curvature 
gravity. To a review of the quadratic curvature gravity see for example the reference 
and the references cited there. 

Consider the quadratic gravitational theory defined by the action 

S = tt-t; ! d^x^{R + + j3R R^^- + IQtiGC^}, (34) 
IottG J 

where Cm is the Lagrangian of matter fields and G the Newton's gravitational constant. 
In writing the action (jH^ we make use of the four dimensional Gauss-Bonnet invariant to 
eliminate the quadratic invariant, R u„n^ R^^''" as is usual in the treatment of the quadratic 
curvature action in four dimensions [15[- 

The field equations derived from 5S = Q are, 

Rf,u - ^9iMuR + aiif,u + PlfMu = SttCT^j,, (35) 

where B.^^, and l^^ are given by 

Hfiu = —2R-^u + '^gfiuOR — -g^uR^ + 2i?i?^jy, (36) 

and 

^^J.l^ = —'2R"^-ua + ORfiu + -^g^iuOR + 2R"^Rau — -gf_iuRaf3R"'^ , (37) 

where □ is the curved space d'Alembert operator. To the theory defined by ()34j) to 
have an acceptable Newtonian limit, the parameters a and f3 must satisfy the following 
constraints |l6|: 

3a + p>0 , p <0. (38) 
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For metrics of the form (^J) we have □ V\, where V^^ is the Laplacian in the transverse 
space. Moreover, by virtue of Q there is no contribution of the Hf^u term and the only 
contribution of the I^y comes from DR^^. Therefore, the field equations reduces to: 

-^[/?Vl + Vl]H{x^, u) = 8nGT^^. (39) 

This equation can be integrated to give 

[/3Vl + l]H{x±,u) = Hi{x±,u) + ah{x±,u), (40) 

where h{x±, u) is a harmonic function of the transverse coordinates, a an arbitrary constant 
and 

VIH^{x^,u) = -IQ'kGT^u. (41) 

Now, we write H{x±,u) as 

H{x±,u) = Hi{x±,u) + H2{x±,u) + ah{x±,u), (42) 

where H2{x±,u) is determined by the equation 

H2{x^,u) = 167iGTuu. (43) 

The functions Hi and H2 are respectively the purely linear and purely quadratic parts of 
the solution. The pp-wave solution to quadratic gravity is given by the metric (Q) in which 
H is given by the equation (jl^ . 

Regarding the equation (|T9|l and substituting in we obtain 

Ao ^ 1^22 = -l^^lH = j^H2{x^,u). (44) 

Thus, the helicity-0 component of the wave is given by the solutions of the equation ()43|) 
multiplied by —1//?, the inverse of the coupling parameter of the R^yR^^ invariant in the 
quadratic gravitational action. We notice that the linear curvature term do not contribute to 
the $22- Therefore, the transverse helicity-0 mode comes only from the quadratic curvature 
terms and is determined by the equation PHj) . This fact contrasts with the result that 
$22 = in empty space pp-wave solutions to Einstein's gravity. There are also also 
the nonvanishing components of helicity-2 namely, A^ and Ay. which can be computed 
from (j2DI), where the effect of the quadratic curvature comes from the if 2 term in (021) • This 
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result indicates that can be a contribution of quadratic curvature terms to the Einsteinian 
polarizations of a null wave. This fact provides a justificative, from a non perturbative point 
of view, to the appearance of a /3 dependent correction in the amplitude of the linearized 
gravitational waves in the transverse traceless gauge when quadratic curvature terms are 
considered as small corrections to the Einstein's general relativity 2]. 



A. An example of impulsive gravitational wave in quadratic gravity 

Let us give an explicit example of a solution obtained in quadratic gravity for a source 
consisting of a finite thin shell of width 2Ro and homogeneous energy density ^o, which 
propagates with the light velocity in the z direction. The T^^j component of the matter 
energy-momentum tensor is: 

Tuu = XQoQiRo-p)5iu), (45) 

where G(x) is the standard step function, p is the radial coordinate from the source axis 
and A is a constant [1]. By solving the equations (jlTjl and with Tuu given by we 
obtain 



Q{Ro-p)S{u) 

0(p-i?o)5(M); 



^j + 4i?o&/i(^ji^o(0 e(p-i?o)5(«)|, (46) 

where k = AttGXqo, K,^ and Ii, are modified Bessel functions and b = \J —fi (see jlj for the 
explanation of the boundary and regularity conditions which are imposed in the derivation 
of the result pUj) ). The cylindrical symmetry of H is due to the symmetry of the source (jlHJ. 
We assume that /3 < since when /3 > there is no acceptable Newtonian limit for the 
nonrelativistic gravitational potential between point masses in quadratic gravity This 
choice excludes the non physical solutions which appear when /? > leading to imaginary 
components in H. Note that the solution is a continuous function of p and diverges 
logarithmically for p ^ oo. Although H is natural quantity to be taken as the amplitude of 
app-wave, this choice is not appropriate to the observational point of view since it contradicts 
the expectation that the wave amplitude must decrease with the distance form the source. 
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Thus, the quantities that can better represent wave amphtudes from the observational point 
of view must be given by the and Aq which determines the relative accelerations 

between the test particles. 

Using the equations (fTUj) and (pUj) we obtain 



Ao = -$22 



1 



e(i?o-p)5H 



(47) 



A, 




tHt 



Ro 
T 



/2(^)e(i?o-p)5(w) 



Q{p-Ro)S{u) 



and 



A. 



1 



Qm'^A = 0. 



(48) 



(49) 



Both Ao and are continuous functions of p and goes to as p — >■ oo. If an observer is 
placed at the source axis (p = 0) with one of his frame direction, for instance the z axis, 
aligned with the propagation direction of the wave, he do not see the helicity-2 component 
in the acceleration pattern of test particles since this observer is located at the symmetry 
axis of the source which has cylindrical symmetry. For this observer y4+ = Ax = and the 
only nonvanishing pattern in the relative accelerations of test particles comes from 



Ao{p = 0) 



2 



Ro ( Ro 



1 



b{u) = -^H{u, 0) = ~H,{0, u), (50) 



which has no contribution from the linear curvature (Einsteinian) part of the theory. If the 
observer keep his frame orientation but, is displaced at a distance p from the source axis, 
he observes only one component of helicity-2, given by PHj) in addition to the helicity-0 one 
given by (jUj). 



V. SUMMARY AND CONCLUSIONS 



We have studied the deviation of geodesies in pp-wave spacetimes by relating the 
nonvanishing NP quantities of a general pp-wave in four spacetime dimensions with the 

13 



Riemann tensor components in a local orthonormal basis. We have showed that the pp-wave 
solutions to quadratic curvature gravity produces relative accelerations between test particles 
located in the geodesies of the spacetime which are transverse to the wave propagation 
direction. These accelerations follow a pattern given by, at most, two components of helicity- 
2, analogous to the polarization pattern of a plane gravitational wave in linearized Einstein's 
gravity, and one of helicity-0. For a general quadratic gravity pp-wave we have obtained that 
there is a helicity-0 pattern in the relative accelerations of test particles which depends only 
on the purely quadratic part of the spacetime metric, namely A particular example of 
an impulsive pp-wave solution to quadratic gravity with cylindrical symmetry was given for 
which we identify one helicity-2 and one helicity-0 nonvanishing components that can be 
observed in the relative accelerations of nearby test particles. The suppression of one of the 
helicity-2 patterns occur due the cylindrical symmetry of the source. For an observer placed 
at p = 0, the helicity-2 component vanishes due to the symmetry of the source. For this 
observer only the helicity-0 pattern, which depends only on the quadratic curvature part of 
the metric, survives. 

The approach by which the relative accelerations of nearby test particles in a local 
orthonormal basis was obtained can be used to obtain the geodesic deviations in more 
general spacetime metrics. An interesting study that can be carried out within the context 
of quadratic curvature gravity concerns the geodesic deviations of test particles in presence 
of non-null wave-like solutions to quadratic gravity such as those obtained in [3]. However, 
we left the investigation of this subject to be carried out in another future work. 
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APPENDIX A 

We transcribe in this appendix the definitions of the NP quantities for a null tetrad basis 
{I, k,m,ffi} and write the Riemann tensor components in the local (observer) basis {ea} 
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in terms of the NP quantities. 

In the NP notation, the null tetrad components of the traceless Ricci tensor (5*^,^ = 
Rfiu — g^ivRl^) can be written in terms of three real and three complex scalars according to 
the definitions P]: 



$00 = \s^,k^k'' = $00 , $01 = \s,uk^m'' = $10 , 
$02 = \s^,m^m' = $20 , $11 = \s^.u{kn'' + m^rff) = $n , 

$12 = \s,,Fm^ = $21 , $22 = ^S^J^r = $22 . (Al) 

The Ricci scalar is denoted by 

^ = = Y^(^(3)(4) - ^(1)(2)), (A2) 

where R{a){b) = Rfiu^(a)^'(b)- "^^^ Weyl tensor components in the null tetrad basis can be 
written in terms of the five complex scalars : 



^0 = C^^^sk^m^k^'m^ , = C^^^sk'^l'^k^m^ 
^2 = \c^,^5k^r{kn' -mm') , ^3 = C^,^d''k''Pm' 



^4 = C^^^sl^'^Pm^ . (A3) 



The Riemann tensor components in the local basis {ea} in terms of the NP quantities 
are: 



^i6i6= 2^^^o+2^^^^~^^^2 + 2$22 + 2*00 - 3fJe$o2-2A , (A4a) 



1 1 

^1620 = 2 ^"^^0 - -5>m^'4 - 5>m$o2 , (A4b) 



= -SRe^i + 3?e^3 - 3?e$oi + 3?e$i2 , (A4c) 



^2026 = -^3f^e^o-^3f^e^4-3fJe^2 + ^$22 + ^$00 + 5Je$o2-2A , (A4d) 
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i?2630 = ~ ^m'^s - '^m^oi + Qm^i2 , (A4e) 



^3030 = 23f^e^2 + 2$ii - 2A . (A4f ) 
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